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Abstract 

We consider the coupling between three dimensional gravity with zero cos- 
mological constant and massive spinning point particles. First, we study the 
classical canonical analysis of the coupled system. Then, we go to the Hamil- 
tonian quantization generalizing loop quantum gravity techniques. We give a 
complete description of the kinematical Hilbert space of the coupled system. 
Finally, we define the physical Hilbert space of the system of self-gravitating 
massive spinning point particles using Rovelli's generalized projection opera- 
tor which can be represented as a sum over spin foam amplitudes. In addition 
we provide an explicit expression of the (physical) distance operator between 
two particles which is defined Dirac observable. 
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1. Motivations 



Loop Quantum Gravity provides a promising proposal of quantum theory of gravity 
PlEj- In contrast to usual quantum field theories, it is a background independent, 
non-peturbative approach to the problem of unification of general relativity and 
quantum mechanics. Loop quantum gravity is based on the canonical formulation 
of general relativity expressed in terms of connection variables. The quantization 
program can be presented in terms of two (technically) distinguishable steps: in 
the first step the constraints governing general covariant dynamics are represented 
as quantum operators in a Hilbert space of kinematical states. The second step 
consists of characterizing the physical Hilbert space of solutions of the quantum 
constraints. The first step is well understood. The kinematical Hilbert space is 
essentially spanned by polymer-like excitations known as spin-networks. Moreover, 
they provide a definition of the notion of quantum geometry of space as they corre- 
spond (in a suitable sense) to eigenstates of geometric operators. The quantization of 
the scalar constraint (governing the quantum dynamics) has been achieved although 
due to certain ambiguities in the regularization several proposals (mathematically 
equally consistent) have been put forward. The ambiguity is expected to be resolved 
by the restrictions that would arise in trying to reproduce the semiclassical limit of 
the theory. In order to achieve this one would need to resolve the second step men- 
tioned above, i.e., to characterize the physical Hilbert space for a given quantization 
of the scalar constraint. This second step remains to a large degree open in loop 
quantum gravity. The spin foam approach provides a systematic tool to describe the 
dynamics aiming at the construction of a mathematically well defined notion of path 
integral for loop quantum gravity tl8|. For a novel proposal to tackle dynamics in 
loop quantum gravity which also leads to a spin foam representation see [^Hl 

Despite the difficulties associated to dynamical issues, important physical ques- 
tions have been successfully addressed by the theory. In particular, loop quantum 
gravity provides a microscopical derivation of black whole entropy ^IT] and predicts 
that space is discrete at the Planck scale [20j. More recently, novel insights about 
the nature of the classical big-bang singularity and the physics of initial conditions 
have been achieved in the context of loop quantum cosmology j3] . The prediction 
of discreteness of space at the Planck scale is certainly the most important result 
of loop quantum gravity; playing a crucial role in the other physical predictions of 
the theory. Although the possibility of direct observation of discreteness seems not 
plausible at this stage, one must ask the question at the conceptual level, i.e., is it 
possible to imagine a physical process that could measure, without any ambiguity, 
geometrical properties of space, as lengths, areas or volumes? Making measure- 
ments of geometrical properties of space-time is meaningless without considering 
the presence of an observer in the theory [221 • I^ practice, this means that one 
has to couple the gravitational field to a physical field associated to the observer in 
order to construct physical geometrical operators or Dirac observables in the Hamil- 
tonian framework. However, this idea is technically very complicated to implement 
in general. Already at the classical level the definition of observable quantities is 
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quite involved. Aside for a few examples that are useful in special cases (e.g. the 
area of black hole isolated horizons [1]) the construction of Dirac observables in the 
quantum theory remains quite open. 

Three dimensional gravity offers a good framework to test these ideas. In the 
absence of external field, it is a topological field theory and therefore only admits a 
finite number of degrees of freedom. This property makes the classical theory exactly 
solvable and different quantization schemes have been proposed in the literature 
[7j. In particular, the loop quantum gravity program can be completely realized in 
that case: the kinematical Hilbert space is well understood and, in a companion 
paper ^H], we make use of spin-foam models to solve explicitely the dynamics. 
Geometrical operators, as length operators, have been studied and their spectrum 
have been computed in both Euclidean [22] and Lorentzian regimes ^UI- These 
operators are not Dirac observables and cannot be, a priori, related to a physical 
process. In the spirit of the previous lines, the simplest way to construct geometrical 
observables is to couple the gravitational field with particles and consider distances 
between these particles. At the classical level, this has been successfully done using 
the Chern-Simons formulation of gravity Distance observables between particles 
are given explicitely, they form a Poisson algebra which naturally exhibits solutions 
of the classical dynamical Yang Baxter equation. This very nice algebraic result 
generalizes the combinatorial formalism of pure Chern-Simons theory [BJ and allows 
to consider quantizing the coupled system by the use of dynamical quantum groups. 

Nevertheless, in this article we will consider the loop quantization of three dimen- 
sional gravity coupled to massive and spinning point particles. First, we construct 
classical distance observables in the Dirac sense. We then proceed to the quanti- 
zation, find their spectrum, and compare it with (non physical) distance operators 
previously considered. We describe the full dynamics of the self-gravitating massive 
spinning particles system in terms of the spin foam representation. Recently, the 
inclusion of particles in the Ponzano-Regge model has been considered in a rather 
different context (ITI . 

The paper is organized as follows. In the next section, we present the classical 
coupling between point particles and the three dimensional gravitational field: we 
adapt the results of [Hj in the case of Euclidean gravity with zero cosmological 
constant and we give an explicit expression for the distance between particles which 
are Dirac observables. In the third section, we quantize the coupled system a la loop 
quantum gravity: we first give a complete description of the kinematical Hilbert 
space, then we construct the physical Hilbert space. To describe the kinematical 
Hilbert space, we extend the notion of spin-network and we define generalized spin- 
network states which include particles. To construct the physical Hilbert space, we 
generalize the construction presented in the companion paper jTH]: we propose a 
regularization for the projector on physical states and we obtain a description of 
the physical scalar product in terms of the spin foam representation. The obtained 
spin foam model generalizes the Ponzano-Regge model in the sense that it includes 
massive and spinning point particles and, as in JE] , it is not defined on a background 
triangulation as all regulators are removed in the definition of the physical inner 
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product. 



2. Classical formulation of the theory 

Leaving out the problem of degenerate metrics, it is well known j2ZI that three 
dimensional gravity is equivalent to a Chern- Simons theory whose gauge group G 
depends on the signature of the metric and on the sign of the cosmo logical constant. 
Furthermore, the coupling between a point particle and the gravitational field can 
be formulated as a minimal gauge coupling [SJ [T3] . 

This section aims at using this formulation to make a precise classical analysis 
of the coupling between particles and Euclidean gravity without cosmological con- 
stant. For that purpose, we will first present a non-usual formulation to describe 
the dynamics of a free Euclidean relativistic particle in such a way that we could 
immediately apply, in a second step, gauge theory methods to describe the coupled 
system. 

2.1. Free relativistic particle 

Our formulation is based on the fact that we can identify the degrees of freedom of 
a massive spinning relativistic particle evolving in the three dimensional Euclidean 
space with an element of the (universal covering of the) isometry group G = 
ISU{2) of (see for example [E]). Indeed, any element X e G can be decomposed, 
by definition, as the semi-direct product X = (A, q) of a rotation A e SU (2) and 
a translation q. Then, we can obviously identify the element q with the position 
of a particle in E^ and we will see latter that A allows to construct a "spin-vector" 
associated to the particle ^3]. In the sequel, we make extensive use of the vectorial 
representation of the group G where A is a 3 x 3 rotation matrix and g is a three- 
vector of E^. Notation and properties of the Euclidean group G = ISU{2) and its 
Lie algebra g = isu{2) are exposed in the appendix. 

2.1.1. Algebraic description 

Given a point particle of mass m and spin s, its degrees of freedom are completely 
determined by an element X = (A, q) of the euclidean group G and its dynamics 
between times ti and t2 is defined by the following first order action [Tl] : 



Sp[X] = dt < x(m, s) , X > where ■§) = + sPq Eg. (1) 



The application <,>: 0X0 —>■ Cisa non-degenerate invariant bilinear form on 
g (appendix). We will introduce the notation C ~ G to denote the configuration 
space of the point particle. In that framework, Euclidean transformations (analog 
of Poincare transformations) of the particle degrees of freedom are defined by the 
following map: 




C X G 



C 



gX. 



(2) 
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Applying these transformations in the vectorial representation of G, one can in 
particular recover the well-known Euclidean transformations of the position q of the 
particle. 

Before going into a precise canonical analysis, let us show in the Lagrangian 
framework that Sp reproduces the motion of a free massive spinning particle on . 
For that purpose, we compute the variation of the action 5Sp with respect to X and 
we see immediately that the variational problem is well defined when the quantity 
< xij^is) , X~~^5X > vanishes on the boundaries t = ti and t = ^2- Actually, this 
condition fixes half of the degrees of freedom of the particle at initial and final times 
and, in the particular case of a purely massive particle [s = 0), it is equivalent to 
the condition 6q = 0, i.e. the position of the particle is fixed on the boundaries. 
With these rather standard boundary conditions the equations of motion are 

[X(m,.), X-i— ] = ^ X-'— = aJo + PPo, (3) 

where [, ] is the commutator in the Lie algebra. 

In the last equation, a and (3 are, a priori, arbitrary functions of time. The 
equivalence (jH} is trivial and general solutions of the equations of motion read X{t) = 
X{ti) exp J^^{aJo+pPo) where X{ti) is an element of G defined by initial conditions. 
In order to clarify the relation between these solutions and the motion of a massive 
spinning particle in E^, let us write them explicitely in the vectorial representation: 

/•t2 

A{t) = A(ti)exp / dtaJo and q{t) = q{ti) + T{t)A{t)Po . (4) 

We have introduced the notation T{t) = J^* pdt. As a result, the dynamics of the 
particle coordinates takes the form = q^{ti) + —p'^T{t) where we have defined 
=< Jj,APo >= ^ < A^^JjA,Po >• Thus we obtain the geodesic equations of 
the metric space E^. Notice that are constants of motion satisfying the constraint 
p2 _ ptpi _ ^2 therefore are interpreted as the components of the 3-momentum 
of the particle. Besides, the function T{t) is viewed as a clock^ associated to the 
particle which measures time evolution. Following the spin vector a of the 
particle is proportional to the momentum, i.e. a* = — and, therefore, its dynamics 
is trivial. Note that the momentum components (and also the spin) are expressed 
in a given frame, called internal frame in the sequel. 

This brief Lagrangian analysis makes a crystal-clear contact between the dy- 
namics induced by the action and the motion of a standard free particle in E^. 
Furthermore, it is immediate to see that the dynamics of a free particle in the three 
dimensional Minkowski space can be formulated in the same way but replacing 
the gauge group G by the Poincare group IS0{2, 1) [T3j. Similar formulation exists 
to describe the motion of a particle in the three dimensional de Sitter space dS^ jSj 

"'^The clock T is an arbitrary function but should be a monotonic function of time parameter 
in order to describe causal trajectories. Actually, it is possible to find a derivative gauge fixing 
compatible with the variational problem and such that T{t) is monotonic (cf 2.1.3). 
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where the gauge group is now replaced by the isometry group of dS^ and we can 
extend this formulation to the Anti-de Sitter case. 

This non-usual formulation will appear soon very interesting to describe the 
coupling with a gravitational field and then to quantize the obtained system. 

2.1.2. Canonical analysis 

This section is devoted to the canonical analysis of the free particle action. Note that 
we will present the analysis in a formulation such that one could easily generalize to 
the Lorentzian, deSitter or anti-deSitter cases for example. The reader not interested 
into the details can have the main ideas of this analysis in the footnote ^at the bottom 
of this page which draws its inspiration from |14j . 

To a general study of the Hamiltonian analysis, it will be convenient first not to 
distinguish the position and spin degrees of freedom of the particle. Therefore, we 
consider a parametrization z = {z^)j=i^...fi G of the group G and we will denote 
by X{z) the element of G parametrized by z. 

We associate to each parameter z^ the canonically conjugated momentum tt/ 
such that the Poisson brackets are given by {ttjjZ'^} = 6j. Using the notations 
introduced in the appendix, it will be convenient to define the matrix P = qj^^ 
where we have considered the variables qi = 7rj|^(a = 0, 2;). In terms of these new 

^The variable X G ISU{2) contains all the degrees of freedom of the particle: the position 
and three angles labeling A G SU{2) which carry the information about the direction of the 
momentum of the particle — as we will show below one of the three angles is a gauge degree of 
freedom while the other two label a point on the sphere of directions of pi . From the properties 
of the Killing form (see appendix), it is immediate to show that the action can be rewritten as 
follows: 

S,[,.A| = £di(,™.f-iTr«A-.f))^ 

In this expression, is a unit vector defined by n'^Ji — AJqA^^ and Tr is the trace in the 
fundamental representation of SU{2). 

This action clearly defines a constrained system whose constraints reads |14j : 

$i = Pi — mrii ~ and = ai ~ srii ~ , 

where pi is the momentum (conjugated to position g*) and ai is the spin- vector whose expression 
in terms of the phase space variables is given in 11. In a sense, one can think of the vector-spin 
components as variables "conjugated" to the angles labeling A. It will be convenient to consider 
the total angular momentum ji = Xi + ai where the orbital momentum is defined hy X = q A p. 
From the set of six constraints, one can extract two first class constraints which read: 

Cp = pip^ — rn? ~ and Cj = pif — ms ~ . 

The first constraint is clearly the usual mass shell condition of the relativistic particle and there- 
fore generates reparametrizations along the trajectory of the particle. The second one generates 
rotations of the frame which conserve the momentum p^ (its orbit corresponds to the third angle 
mentioned above). It can be also viewed as a definition of the spin s of the particle. The analy- 
sis of the second class constraints leads to the Dirac bracket whose expression is given in |14| or 
equivalently by the expression (|10() in the core of the text. 
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variables, the Poisson brackets read: 



{Xi,X2} = 0, {P^,X2} = -h2X2 and {Pi, P2} = [^12, ^i] • (5) 

In these equations, we have used the standard tensorial notation and t G g*^^ denotes 
the Casimir tensor associated the the bihnear form <, >. 

The free particle action clearly defines a constrained system and the con- 
straints analysis is done along the same lines as in the de Sitter case In par- 
ticular, conjugated momenta are fixed in terms of the coordinates by a set of six 
primary constraints which can be expressed in the following matricial form: 

C[X,P] = X-^PX + x{m,s) ~ 0. (6) 

Once the primary constraints identified, we introduce the Lagrange multiplier /i G 
C°°(M, 0) in order to define the total Hamiltonian by Htot =< fJ',C[X,P] > for the 
canonical Hamiltonian of the theory is identically zero. Using the immediate identity 
71 1 = — < P, MM~^ >, we see that the free particle action ((H) is equivalent to the 
following first order action: 



Stot[X,P,fi] = - dt P,^M-' > ~Htot[X,P,fi]^ 



(7) 



The total Hamiltonian Htot is necessary to define temporal evolution of the dy- 
namical variables and we show that conservation of primary constraints under time 
evolution does not introduce new constraints. However, the Lagrange multiplier fi is 
restricted to be an element of the Cartan subalgebra, i.e. fi = fipPo + /ijJo- As a re- 
sult, among the six constraints ©, Cp =< Pq, C[X, P] > and Cj =< Jo, C[X, P] > 
are first class. On the other hand, the functions Ci =< C[X, P] > where 
^/ G {Pi, Ji)j=i,2 form a set of second class constraints whose Dirac matrix is given 
by: 

{Ci,Cj} = Ajj =<x(m,5),[6,0] > • (8) 

We have introduced the notations m = m + Cp and s = s + Cj. Note that it is 
crucial to invert the Dirac matrix strongly in order to obtain a Dirac bracket which 
satisfies, in particular, the Jacobi identity. This is the reason why we don't solve 
explicitely the first class constraints in the expression (jH)) of the Dirac matrix. Taking 
this remark into account, we compute the inverse Dirac matrix and we obtain: 

(A^^)^-' =< K{m,S),[^',^^]> where K{m,s) = ""^^^^2^° • (9) 

As a consequence, it is possible to compute Dirac bracket between any two functions 
ip and ip defined on the phase space. In particular, Poisson brackets between matrix 
elements of P and any function on the phase space are not modified by the Dirac 
procedure for P commutes with the second class constraints. On the other hand. 
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Poisson brackets between matrix elements of X are replaced by the following Dirac 
bracket: 



{Xi,X2}d = X,X2ru{m,s) . (10) 

As a consequence, the Dirac bracket appears to be quadratic and introduces the 
r-matrix r : — {(0,0)} ^ Cg> g defined, for any non zero couple (m, s), by 
the expression ri2(m, s) = [K2{fn, s),ti2]. The antisymmetry property of the Dirac 
bracket implies that the symmetric part of r vanishes, i.e. ri2 + = 0. Moreover, 
the Jacobi identity is equivalent to the fact that r is a solution of the classical 
dynamical Yang-Baxter (CDYB) equation associated to the Lie-algebra fl E]- 
This result is particularly interesting and underlines an unexpected link between the 
dynamics of a free relativistic particle and solutions of the CDYB equation. This 
relation could be the starting point of a quantization program for the relativistic 
particle involving quantum groups but we will adopt a more conservative point of 
view in this article. 

From the Dirac bracket (fTUI) . one can establish the equations of motion of the 
dynamical variables in the Hamiltonian framework and one obtains: 

dX dP 

— = {X,Ht,t]D = X/i and — = {P,//^^}^ = 0. (11) 

As one could expect for a free particle, the conjugated momentum P is a constant 
of motion. If one decomposes it as follows P = pVj + j*Pi, then (p*) and (j*) are 
respectively the momentum and the total angular momentum of the particle whose 
expressions are explicitely given by the formulae: 

pVi = -mA Jo A"^ and fPi = +mAJoA"^ g - sA Pq . (12) 

As we could expect, the total angular momentum is the sum of the angular momen- 
tum X = q A p and the spin-vector a = ^p. Using the two invariant bilinear forms 
on Q (appendix), one can determine the relativistic invariants of the free particle as 
follows: 

<P,P>= 2p-j = 2ms and (P, P) = p-p = . (13) 

The dot ■ denotes the usual scalar product on the three dimensional Euclidean space 
E^. As for the equation of motion for X ()11|). it is the same as the one obtained 
in the Lagrangian framework (jH)) when one replaces a and (3 respectively by /xj and 
lip. As a result, general solutions of the equations of motion read: 

X{t) = X(ti)exp [ dufi{u) and P(t) = X(ti)x(m, s)X(ti)-^ (14) 
Jti 

In the vectorial representation, it is clear that these solutions reflect the dynamics 
of a free relativistic particle in E^. At this level, fi appears to be any function of the 
time with values in the Cartan subalgebra of q and thus there is no restriction on the 
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evolution of the particle position along its straightline trajectory. For instance, it is 
possible to find fi such that the particle oscillates along its trajectory and then its 
motion violates causality. In fact, the freedom to choose any particular evolution for 
the particle comes from the invariance of the system under time reparametrizations. 
As a consequence, one can fix this problem by an appropriate gauge fixing 121] ■ 
For that purpose, it is necessary to make a study of the symmetries of the system. 

2.1.3. Symmetries and time reparametrization 

Among different types of symmetries, we distinguish gauge symmetries (local) and 
Noether symmetries (global). In the Hamiltonian formulation, the former are gen- 
erated by first class constraints whereas constants of motion generate the latter. 

In our formulation of the relativistic particle, the generator of infinitesimal gauge 
transformations is given by G[e] =< e, C[X, P] > where e is an infinitesimal element 
of the Cartan subalgebra of q and gauge symmetries on the dynamical variables read: 



When we evaluate the previous Poisson brackets, we assume that the infinitesimal 
parameter e does not depend on the dynamical variables. However, it appears that 
the action ((T]) transforms as 6eS =< ^(^2) — £(ti),x{™j > under these transfor- 
mations and therefore is generically not let invariant. As a result, the invariance of 
the action is obtained if one imposes the vanishing of < e, xi^^ s) > on the bound- 
aries. In that case, the total action (|7j) is invariant under the transformations ()15|) 
if and only if the Lagrange multiplier /i transforms as 6sfi = —e. At this point, it is 
very interesting to note that, if one defines the function M G C°°{[ti,t2],G) by the 
expression 



then the quantity X{t)M{t) is invariant under gauge symmetries and therefore is an 
observable of the system. In fact, from solutions of the equations of motion ()14|). we 
see immediately that X{t)M{t) = X{ti) is fixed by initial conditions. In the sequel, 
M{t) will be referred as the dynamical function of the particle for it contains all the 
dynamical information. 

To clarify the relation between these gauge transformations and usual symme- 
tries of the free relativistic particle, let us remark that SeX = —XX — X{Xfi — X) 
when e = Xfi and A is an arbitrary function of t. As a result, one first notes that 
gauge symmetries are equivalent, up to the equations of motion, to reparametrization 
symmetries of the free particle action (jT)), i.e. t ^ t — X{t). However, gauge trans- 
formations contain more than the usual reparametrization symmetry of the particle 
trajectory for they are generated by two different first class constraints which can, 
in fact, be rewritten as Cp = —p^ —m and Cj = —p- j — s. Thus, one can show that 
the former explicitely implements time reparametrizations of the particle position q 



S,X = {G[e],X} 



Xe and = {G[e],P} 



. 



(15) 




(16) 
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whereas the latter generates rotations of the internal frame which let invariant the 
momentum p. 

We can ask the question of a gauge fixing compatible with the boundary condition 
satisfying by the parameter e. The derivative gauge fixation < /i, xi^^ ^) >= is a 
good candidate. To be more precise, this is a partial gauge fixing and the remaining 
gauge freedom can be fixed using a canonical gauge. In the purely massive case, it 
implies the fixation of the component /ip of the Lagrange multiplier and therefore 
the clock T{t) associated to the particle is a constant (so monotonic) function. 

Before going to the coupled system, it is important to note that the free particle 
action ([Q) is invariant under the global transformations (j21). This symmetry reflects 
the freedom to choose the frame (origin and axes) of the Euclidean space E^. The 
generators of these infinitesimal symmetries are given by the constants of motion 
p\j^ and we have: 

{p\X}d = -m and {J\X}D = -J^X. (17) 

Moreover, and j* commute with the constraints and therefore are observables of 
the theory. It is easy to show that their Poisson algebra satisfies the Euclidean 
algebra, i.e.: 

{p\p'}d = , {f^j'jn = and {p\f}n = ^\v' ■ (18) 

The coupling to the gravitational field consists precisely of gauging the global sym- 
metries fll7|l using, as usual, a G-connection A. However, a naive minimal gauge 
coupling introduces ambiguities and therefore it is necessary to consider a regular- 
ization. This point was completely understood in the de Sitter case and next 
section is devoted to recall the basic ideas in our case without entering into the 
details. 

2.2. Coupling to the gravitational field 

The fact that there is no gravitational interaction between point particles in the 
three dimensional case allows one to study the coupling to each particle with the 
gravitational field independently. For this reason, we will endeavor to present the 
case of one particle coupled to gravity; the generalization to an arbitrary number 
of particles will be immediate. We will consider a three dimensional manifold M. 
whose topology is given by = S x [^1,^2] and we will denote by the trajectory 
of the particle. 

2.2.1. Minimal gauge coupling and regularization 

Up to boundary terms, the first order action of three dimensional Euclidean gravity 
defined on the manifold M. can be formulated as a Chern-Simons action whose gauge 
group is G = 1SV{2): 

Scs\A\ = e^'p [< A^, d,A, > +^ < A^, [A, Ap\ >^ . (19) 
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The connection A = A^dx^^ is related to the spin connection A* = A^^dx^ and the 
triad e* = e^dx^ by the expression A^ = e^^Pi + A^^Ji. This action defines a gauge 
theory and infinitesimal gauge transformations of the connection are given by: 

V e G , s) , 6^A^ = V^i = d,^ + [A„ e] . (20) 

We can use the Chern-Simons connection to gauge the global symmetries of the free 
particle (fTTj) such that the coupled action reads: 

S,\A,X] = dt < x{m,s),X-'— + X-'AtX > . (21) 

The integral is defined along the trajectory of the particle and is, by construction, 
invariant under local transformations (fTTj) when the connection transforms as (|2()j) . 
Finally, the dynamics of the whole coupled system (particle and gravity) is governed 
by the action: 

S[A,X] = Scs[A] + S,[A,X] . (22) 

Equations of motion are obtained by varying this action with respect to the dynam- 
ical variables. In particular, the time component of the connection At appears to 
be a Lagrange multiplier which imposes the following constraints on the curvature 
of the Chern-Simons connection A: 

V . e E. = p s(. - ^ { : ^: ^ -_ (23) 

When we decompose l.h.s. equations (j^Hj) into its translational and rotational parts, 
we obtain well-known constraints (r.h.s.) on the torsion tensor T[e,A] and the 
curvature tensor F[A]. As we could expect, we find that the connection is singular at 
the location of the particle: this is the signature of the well-known conical singularity 
induced by a particle in three dimensional gravity 0. However, the equations 
can only fix the conjugacy class of the holonomy of the connection along a loop 
surrounding the particle unless we specify a given origin on the surface S. This 
equation is then ill-defined and it is necessary to reformulate the coupling (j2H) in 
order to regularize it. Different schemes already exist El and we will adopt the 
one developed in the article that the reader is invited to read if she is interested 
by the details. 

The basic idea consists of replacing the point particle by a small loop £ whose 
time evolution draws a cylinder B = i x [ti,t2] in the space-time A4. In such a 
picture, the coupling is defined as an integral over B which trivially generalizes the 
action (j2ip and we have: 

SrnA,X] = ^ [ dtdip<xim,s),X-'^+X-UtX> . (24) 
Ztt Jq dt 

Note that G [0, 27r[ is the angular variable parametrizing the loop £ and the origin 
<y9 = is arbitrary. In this action, X still depends only on t and we see that the 
coupling is, in fact, defined with the mean value of At along the loop £. 
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The dynamics of the gauge field A is, a priori, described by a Chern-Simons action 
(fT^ . However, because of the regularization, the space-time admits a boundary B 
which breaks the gauge invariance. To restore the invariance, first we have to ask 
the question of what type of symmetry we would like at the boundary. We will 
naturally extend the gauge transformations (j2Uj) to the case where the manifold 
Ai admits a boundary but we will require that the infinitesimal gauge parameter ^ 
evaluated on any point of B does not depend on the angle ip which does not have any 
physical meaning. Finally, to have a theory invariant under these transformations, 
the Chern-Simons action have to be supplemented with the following boundary term: 

Sb[-A] = / dt dip < At, > . (25) 
Jb 

It is trivial to verify that the resulting action ^csf^] + Ss[A] is invariant under the 
required gauge transformations. Then, the regularized action for the coupled system 
is given by: 

S'^'^iAX] = Scs[A] + Ss[A\ + S:'^[A,X] . (26) 

The generalization to an arbitrary number of particles (pj)i=i,... ,Ar is straightfor- 
ward: one associates a boundary Bi = ii x [ti,t2] to any particle whose dynamics 
is governed by the action ()24|) and one adds as many boundary terms (P^jl as there 
are particles to fix the problem of the symmetry breaking. The dynamical variable 
of the particle pi of mass and spin Si will be denoted by Xj = (Aj, g^) G G. 

2.2.2. Hamiltonian analysis: symmetries and observables 

This section aims to present the main results of the canonical analysis of the coupled 
action (j^ . The reader interested in the technical aspects is invited to read the 
article In particular, the constraint is regularized and is replaced by the 
following constraint defined for any v G C°°(S, g): 

$(t;) = / d^x e"^ < V, Tab > +2 [ d^ <v,A^ + ^Xx{m, s)X-^ > . (27) 
is Je 27r 

As a consequence, the connection is fiat on the surface S but its holonomy Hg 
around the loop i is given by = Xexp(— x(m, s))X~^. Note that there is no 
more ambiguity and the singularity has been removed: the connection A is now 
defined everywhere on S. 

When V is constant along i, then the constraint ()27|1 is first class and therefore 
generates gauge transformations. Besides, the canonical analysis shows that the 
constraints (jH)) remain and G{e) =< e,C > are still first class constraints when 
e is an element of the Cartan subalgebra of g. Let us denote by a = {e, v) the 
set of gauge transformations parameters and the gauge transformations generator 
is defined by T{a) = G{e) + ^(f). Thus, the symmetries generated by r(a) on the 
variables X and on the spatial components of the connection {Aa)a=i,2 read: 

6aX = -Xe-v{i)X and d^Aa = VaV (28) 
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Figure 1: Example of a generalized spin- network observable. 

where v{i) denotes the evaluation of the function v on the loop i. These results are 
trivially extended to the case where there is an arbitrary number of particles and 
we will denote by the infinitesimal parameter of the reparametrization symmetry 
of the particle pi and v{£i) the evaluation of v on the loop £j. 

In the absence of particles, the theory reduces to a Chern-Simons theory asso- 
ciated to the group G and the physical phase space is known as the moduli space 
of flat G-connections. Classical observables are defined in terms of spin-networks 
whose Poisson algebra is well-known ^S]. Different attempts to quantize the theory 
has been developed the last fifteen years but the combinatorial quantization method 
is indisputably the most powerful (see [Oj and references therein). This quantization 
scheme uses quantum groups as a cental tool. 

The combinatorial description of Chern-Simons theory can be generalized in the 
presence of point particles [S]. In particular, we can define new observables which 
capture the dynamics of the particle and contain, in particular, the relative positions 
between particles. To understand this point, let us define the following function: 

0^[A,X,,X^] = M,{t)-'X,{t)-%[A]X,{t)Mj{t) . (29) 

In this definition, [/^[^] denotes the holonomy of the connection along the curve 7 
whose starting point and end point are respectively the points of coordinates ipi G ii 
and (fj G ij (figure [T]); Mj G C°°{[ti, ^2], G) is the dynamical function of the particle 
pi (fTB|) . As the connection A is fiat, depends only on the homotopy class of the 
path 7. Moreover, it is clearly invariant under the gauge transformations and 
therefore is an observable. Similarly to the free particle case, equations of motion 
imply that Xi{t)Mi{t) = Xi{ti) represents the initial configuration of the particle pi: 
this is the reason of the invariance of the observable (j29p under reparametrizations 
of the particles. 

The observable 0^[A, Xi, Xj] can be interpreted as the (initial) configuration 
variable of the particle pj measured in the rest frame of the particle pi. To be more 
explicit, we shall write 0^[A, Xi, Xj] in the vectorial representation and we obtain: 

/ A,(ti)-i/i^[A]A,(ti) A.,{h)-\h^[A]qj{ti) + q^e, A] - qi{h)) 

V 1 

h^[A] = Pexp A denotes the holonomy of the spin-connection A along the path 7 
and we have introduced the notation g^[e, A] = f^h^^x[-A]ea{x)dx°' where the path 
J < X (resp. 7 > a;) is the part of 7 which ends to (resp. starts from) the point 
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X G 7- Therefore, the translational part q-y{i,j) of (jHUI) clearly measures the 
position of the particle pj in the rest frame of pi Then, the relative distance 
T>^{i,j) is given by the formula: 

'^'^{hj) = q-fihjY q-f{ij) ■ (31) 

The rotational part A^{i,j) of represents the (rotation) transformation which 
sends any observer in the rest frame of the first particle to any observer in the 
rest frame of the second one. In particular, it allows to compute the momentum 
components p^(i, j) of the particle pj in the rest frame of pi as follows: 

p^(2,j) = rrij < Pk , A^{i,j)JoA^{t,j)-^ > . (32) 

We will denote by Prest(0 ^^e components of the momentum of pi in its rest frame. 

The distance operator (jHTj) is a gauge invariant function which depends only 
on the homotopy class of 7. Furthermore, when the particles are purely massive 
(sj = 0), V^{i,j) does not depend on the choice of the starting point ipi G ii and 
the end point ipj G ij. When particles possess non trivial spins, the local geometry 
around one given particle pi says that two different points of ii are identified but 
at different time. As a result, one does not expect the distance T>^{i,j) to be 
independent of ipi and (pj, which is indeed the case. 

Let us remark that we deal with the particles reparametrization invariance by 
considering initial configuration variables in the definition of the observables ()29|) 
which makes use of the Lagrange multipliers /ij. Therefore, these observables do 
not capture the dynamics of each particle. Actually, it is necessary to consider a 
gauge fixing to recover the dynamics in the same way as we did in the free particle 
case (cf 2.1.3). However, it is possible to construct non trivial observables without 
considering any choice of initial conditions. A first example of such observables is 
given by the following matrix element: 

e^{t,j) = <k I Ai{t)-^h^[A]A2it) I eo> . (33) 

Notations are recalled in the appendix. It is straightforward to show the identity 
9^{i,j) = p^{i,j) ■ Prestii) and thus (jH^ is physically interpreted as the Euclidean 
boost parameter between the two particles. An other interesting example reads as 
follows: 

d,{^,J) = '^i:'^"V^f'^"^^^-f^^' . (34) 

II P^[hj) /\Prest{l) II 

This expression gives the minimal distance between the trajectories of the two parti- 
cles, is clearly gauge invariant and does not depend on the choice of an initial space 
shce. 

The analysis of second class constraints is done along the same lines as in jS]. 
One can compute the Dirac matrix, inverse it and therefore get the Dirac bracket. 
It appears that the Dirac bracket between spin-networks which do not touch the 
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boundaries is still given by the Goldman symplectic structure [O]. Moreover, the 
Dirac bracket between particles degrees of freedom is not modified and is still given 
by (fTUI) . However, the Dirac bracket between functions depending on the boundaries 
connections is generically very complicated. In the sequel, we will refer to the Dirac 
phase space Vd the set of classical states once the secondary constraints have been 
implemented. Any point x G Vd is completely determined by a G-connection A and 
a sequence (Xi, ■ ■ ■ , X^) of elements of G. 

The set of observables ()29|) forms a Poisson algebra whose structure is very 
nice and introduces solutions of the CDYB equation jH]. This result motivates a 
study of the quantization using combinatorial formalism techniques. However, the 
method is far from being straightforward and needs to be developed to understand 
first representation theory of non-compact dynamical quantum groups and then to 
construct unitary irreducible representations of the observables quantum algebra. 
This program clearly deserves a complete study. 

The following aims to explore loop quantum gravity techniques as an interesting 
alternative of the quantization. In particular, we define the notion of kinemati- 
cal Hilbert space of the coupled system, we quantize the distance observables and 
compute their action on particular states. Then, we propose a spin-foam picture to 
define the physical scalar product generalizing the results obtained in the companion 
paper [TB] . 

3. Loop quantization 

Loop quantum gravity is an attempt to describe a background independent quanti- 
zation of general relativity (see PP for a review). This program clearly distinguishes 
two different steps: the description of the kinematical Hilbert space in terms of 
spin- networks and the analysis of the dynamics using spin-foam models |18j . 

This section aims to adapt loop quantization techniques in the case of three 
dimensional gravity coupled to point particles. For that purpose, we decompose 
the first class constraints $ (j27|) into its torsion part $7- and its curvature part $0 
defined by: 

^t{v) = ^{v'Ji) and ^c{v) = ^{v'Pi) (35) 

for any element v G C°°(S,M'^) which is constant on the boundary i. These con- 
straints respectively generate the rotational and the translational gauge transforma- 
tions (PHj) on the dynamical variables. In fact, they are regularized versions of the 
ill-defined constraints (r.h.s. 1221) • 

The strategy will consist of quantizing the theory before implementing the first 
class constraints. The kinematical Hilbert space TYq will be defined as the set of states 
satisfying all first class constraints but not the curvature constraint $c endowed with 
a suitable scalar product. Then, we will define the physical Hilbert space Hphy from 
the kinematical one by introducing a projector which implements the constraint $c- 
Finally, we will make use of this projector to define a spin-foam model describing 
dynamics of massive spinning particles coupled to the gravitational field. 
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3.1. Kinematical Hilbert space 

The first step consists of quantizing tlie Dirac pliase space Vd endowed witli its 
Dirac symplectic structure. Ttie Dirac bracket between particles degrees of freedom 
is given by the quadratic bracket (jl(J|) and the Dirac bracket between geometrical 
degrees of freedom (not leaving on the boundaries) can be written as follows: 

{e^(x) , A[{y)}n = e,, r^^' 5{x ~ y) , W x,y e S - • (36) 

All the other Dirac brackets vanish but the Dirac bracket between geometrical de- 
grees of freedom leaving on the boundaries of S which are generically very com- 
plicated. However, we won't consider, in this article, this boundary Dirac bracket 
which, nevertheless, deserves a detailed study we postpone for future investigations 

nzi. 

As a result, we can quantize independently geometrical and particles degrees 
of freedom. Loop quantum gravity provides a straightforward quantization of the 
geometrical part (jSHI)- The quantization of the particle degrees of freedom is pre- 
sented in the first subsection. Then, we describe the kinematical Hilbert space of the 
coupled system. Finally, we compute the action of distance operators on particular 
kinematical states. 

3.1.1. Quantum relativistic particle 

At the kinematical level, the particle p is classically characterized by its position 
q and an element A G SU{2). The system is quantized a la Schrodinger and we 
consider the polarization such that A is the configuration variable. This choice of 
polarization is the analog of the momentum polarization in the usual formulation of 
the free relativistic particle. 

We will denote by Pol [A] the vector space of polynomial functions of A whose a 
basis is given by matrix elements of unitary representations of SU{2), i.e.: 

Pol[A] = ©igij^Pol,[A] where Pol, [A] = {n{A)l \ l,n e . (37) 

Notations and basic notions on representation theory of SU{2) are presented in the 
appendix. We endow this vector space with a Hilbert structure making use of the 
SU{2) Haar measure dfi as follows: 



Wf,g e Pol[A] , <f,g>^ j dfi{A) f{A)g{A) . (38) 

We will denote by Ho,i{p) the Hilbert space obtained by completing PolJA] with 
respect to (j38|l and the resulting Hilbert space 7io(p) = ©t'^o,j(p) will be called the 
free particle kinematical Hilbert space in the sequel. 

At this point, we can define an action of quantum observables on ?-^o,i(p)- For 
that purpose, let us recall that the quantum algebra of observables is generated 
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by the elements p'^ and j'^, k = 0,1,2, which satisfy the following commutation 
relations: 

iP" , P'] = , , ]'] = e'^nPn , []' , j1 = e'^nJn • (39) 

From the choice of the polarization, it is natural to think that p'' acts by multipli- 
cation and j'' acts as a derivation in the free particle kinematical Hilbert space. To 
be more precise, one can check that the following action on Ho,i{p)- 

p'> ^(A)l = / HaY^ = -m<Po,A-V,A>^(A)t,, (40) 
r > = ^{A)l = niU , (41) 

defines a unitary irreducible representation of (jHH)- We have introduced the notation 
for the left-invariant derivation on the group SU{2). For any i e ^N, the basis 
(jTTj) of Pol, [A] diagonalizes the angular momentum whose eigenvalues are obtained 
from the following actions: 

jo> ^(A)i = / ^(A)^ and f> k^t = ^i^ + I) k^t ■ (42) 

To define the physical Hilbert space of the system, a preliminary work consists 
of quantizing the constraints and computing their action on ?io(p)- The constraint 
Cp = -^p^ — m is quantize without ambiguity for the components of the momentum 
commute. Thus, we have Cp = ^p^ — m and it is straightforward to see that: 

V/g7^o(p) , C'p> / = 0. (43) 

As for the constraint Cj = ^p ■ j + s, it is also quantize without any ambiguity and 
a direct calculation shows that the action of Cj = ^p ■ j — s on Pol [A] reads: 

Vze^N and V/,nGH,i] , Cj\>^{K)l = {n ~ s) ^(A)^ . (44) 

As a result, the vector space of physical states is the sub-space of Pol [A] (jH7|) defined 
by the following direct sum: 

Pol[Ar ^ 0{^(A)^|/ e H,^]}. (45) 

i-s eN 

As we could expect (in the Euclidean case), we remark immediately that the spin 
of the particle s must be half-integer. Moreover, the direct sum over SU (2) finite 
dimensional representations is restricted to representations i such that i — s 
is a positive integer. In fact, this result is a straightforward consequence of the 
fact that the total angular momentum j is the sum (fT^ of the vector spin a and 
the orbital momentum A which admits integer eigenvalues in the quantum theory. 
In the case of a spinless particle (s = 0), the vector space of physical states 
can be trivially identified to the space of polynomial functions on the two-sphere 
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5*^ = SU{2)/U{1), i.e. the space of spherical harmonic functions. The SU{2) Haar 
measure naturally endows (jl^ with a Hilbert structure and the Hilbert space H{p) 
obtained by completion is precisely the physical Hilbert space of the free relativistic 
particle. 

The quantization of a system of an arbitrary number N of non-interacting rel- 
ativistic particles V = {pk\k = 1, ■ ■ ■ , A^} is straightforward. The physical Hilbert 
space of the system H(V) consists of the tensor product of the physical Hilbert space 
of each particle, i.e. H(V) = ®kH{ipk)- 

Before going to the coupled system, it is interesting to note that the previous 
description gives a "background independent" description of the dynamics of rela- 
tivistic particles in the sense that we do not really make use of the metric structure 
of E"^. Given a manifold = S x [ti,t2]; the dynamics of each particle pk is char- 
acterized by a worldline drawn on M.. In the Hamiltonian framework, each particle 
pk is associated to a point G S endowed with a mass ruk and a spin Sk- As a 
result, a physical state can be described by a set of points on E colored by a couple 
{ik,h) £ |N X [—ik,ik]- Such a description is clearly the analog of spin-network 
states for quantum gravity. 

At this point, one can ask the question whether it is possible to extract geometric 
informations (as distance between particles) from such a "spin-network state". A 
priori, the distance between two particles pk and pn is clearly not an observable 
and therefore its action on the physical Hilbert space T-C{V) is not well-defined. 
However, the minimal distance dmin{k, n) between two worldlines is a gauge invariant 
quantity (i.e. invariant under reparametrizations) and obviously contains geometric 
informations; its expression reads: 

dmin{k,n) = ^ \{qk - qn) ■ {Pk A Pn)\ = ^ \jk ' Pn + jn ' Pk\ ■ (46) 

It is straightforward to quantize this observable, to compute its action on the two- 
particle state vr {AkY^ (8> (An)s^ and to notice that it is generically not diagonal. 
In the case where the particle pn is in rest, i.e. A„ = 1, then the action simplifies 
drastically and is given by: 

dmin{k,n) > 'vr(Afc)t'; ® I = h {Ak)% ® 1 • (47) 

The choice A^ = 1 physically corresponds to consider the particle pn as a classical 
observer who makes measurements on the quantum system. In particular, its mini- 
mal distance to the particle pk is quantized and is given by the magnetic number Ik 
of the particle. 

3.1.2. Quantum kinematical states 

This subsection aims to give a construction of the kinematical Hilbert space of the 
coupled system in analogy with the pure gravity case [T^ . 

First, let us recall the definition of the space of cylindrical functions Cyl(S) on 
the surface E. Given on oriented graph F immersed in E, we denote the set of its 
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edges and its vertices respectively by Sr = {71, ■ ■ ■ ,7£;r} and Vr = {fi, ■ ■ ■ ,vv^}- 
An element "^rj G Cyl(S) is a complex valued functional of the (generalized) 
spacelike connection A labeled by a finite graph F G S and a continuous function 
/ : 5f/(2)>^^r ^ c and defined by: 

^rj[A] = fih,,[A],-.. ,h,,^[A]). (48) 

We recall that hy[A] denotes the holonomy of the connection along the path 7. Given 
two cylindrical functions ^E'riji and ^E'r2,/2 and the minimal graph T containing Fi 
and F2 as subgraphs, then one can trivially extend these two cylindrical functions 
to cylindrical functions with support on T fT^; we will respectively denote by /12 
and /21 the trivial extensions of /i and /2 on T. As a result, one can easily endow 
Cyl(S) with a Hilbert structure using the SU{2) Haar measure d/i as follows: 

/Er 
n Mh^) fi2{hi,--- ,hEr) f2i{hi, ■■■ , hEr) ■ (49) 
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The completion of Cyl(S) under this inner product defines the auxiliary Hilbert 
space Tiauxi^)- This Hilbert space provides a quantization of the classical Poisson 
algebra ()36|) in the case where E has no boundary. To be more precise, it is a 
quantization of the algebra of holonomies of the connection along any loop on S and 
the triad. The former acts by multiplication and the latter as a derivative operator 
whose action is formally given by the following differential operator: 

Our concern here is not to enter into the details of the definition of its action on 
any cylindrical function which has been extensively studied in the literature (see JUI 
and references therein). Nevertheless, it will be useful for our purpose to recall its 
action on a given holonomy h^[A] which reads: 

e^(x) x[C]" > h^[A] = s{x,C,l) h^<x[A] Jk h^>M] ■ (51) 

In this expression, we have denoted x[Q = x[(]°'da the tangent vector to a curve ( 
at X] the quantity s{x, C, 7) = if x is not a point of 7; otherwise it is given by the 
determinant s{x, (, 7) = ea6x[C]"a;[7]^. As a result, it takes value on the set {—1, 0, 1} 
depending on the relative orientation of the curves 7 and ( at the point x. 

The kinematical Hilbert space of pure gravity TCo(J2) is the subspace of gauge 
invariant functions of HauxiX')- Gauge invariant spin- network states provide TCo{T,) 
with an orthonormal basis. The reader is invited to look at the review paper pP to 
have a precise derivation of the kinematical Hilbert space in the four dimensional 
case. 

The next step consists of constructing the kinematical Hilbert space of the cou- 
pled system. For that purpose, we generalize the notion of space of cylindrical 
functions to the following vector space: 

Cyl(S, V) = Cyl(S) ® Pol[Ai]'"" ® ■ ■ • ® Pol[Aiv]'™ • (52) 
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Figure 2: Example of a boundary vertex. The boundary i is associated to a representation 
j£ virtually represented by an extra link oriented toward the boundary; the magnetic 
number associated to the departure point is fixed to the spin s of the particle. The 
boundary is also colored with an intertwiner li from in-coming representations (containing 
the virtual one) to out-going representations. 

The inner products on Cyl(S) and on Pol[A] ()38|) naturally endow the space of 
generalized cylindrical functions Cyl(S,'P) with a Hilbert structure. The resulting 
Hilbert space Hauxi^^'P) will be called the generalized auxiliary Hilbert space. 

The generahzed Gauss constraint dSHl) defines SU{2) gauge transformations 
on the variables of the coupled system. In particular, finite gauge transformations 
of the holonomy h^[A\ of the spin-connection A along an oriented path 7 and of the 
matrices Aj e SU{2) read: 

h^[u] I — > g{s^)h^[A]g{t^Y^ and Aj 1 — > g{t^l^i . (53) 

In these expressions, g G C°°(S, 5^7(2)) whose values g{li) on each circle £j is con- 
stant; and are respectively the source and target points of 7. 

By duality, the action (jSHj) induces an action (a co-action to be more precise) 
on Haux(X','P) and the kinematical Hilbert space of the coupled system Ho(T,,V) is 
defined by the Hilbert subspace of gauge invariant functions of Hauxi^yV). 

In order to provide 7io(S,P) with an orthonormal basis, we need to introduce 
the notion of generalized spin-network. For that purpose, we consider an oriented 
graph r immersed in the surface with boundaries E. Each edge 7 G is colored as 
usual by assigning to it a finite dimensional irreducible representation of SU{2). 
Any vertex f G Vr such that v does not belong to the boundaries Ui£i is associated to 
an S'[/(2)-intertwiner from in-coming representations to out-going representations. 
So far, everything is similar to the usual definition of spin-network but the novelty is 
that the boundaries of S are also colored. Indeed, we associate to each boundary £j a 
couple {ji^, t^J of a finite dimensional irreducible representation j^. and a intertwiner 
Le- of SU{2) in the way presented in the picture 

A generalized spin-network is defined by = {F, {i'y)'y^eTi {'-v)v&Vti Uin '-^J*}- To 
each generalized spin-network S, we associate an element G Ho(T,,V) (called a 
generalized spin-network state) as follows: 

N N . 

^s[A;{AiU ^ ((g)^.J((g)^.)((g)'^'(A,,)J((g) 'v?(^[A])) . (54) 
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Figure 3: The two particles spin- network state. Note that the intertwiners associated to 
the two particles v^^ and vi^ are different. 

This definition is a straightforward generahzation of the usual one. Note that 
there is one free magnetic number associated to each particle link 7^. which contains 
the information of the spin Sj of the particle pi. It is clear that \1/5[A; (Aj)j] is 
invariant under S'f/(2)-gauge transformations (j53|l and therefore is an element of 
T-Cq{T,,V). Moreover, the set of generalized spin-network states form a complete 
orthonormal basis in the generalized kinematical Hilbert space Ti.o{^, V). This states 
coincide with those defined in JT] on a fix triangulation. 

As an example of generalized spin-network state involving particles degrees of 
freedom, let us write the two-particles state \l/^j[y4; Ai, A2] defined by (figure Ej): 

^^,,-[A;Ai,A2] = <e''\n{A^'h.,[A]A2)\es,> . (55) 

In this formula, 7 is an oriented path linking the particles pi and p2- Physically, 
this kinematical state represents a system of two self-gravitating particles with the 
same total angular momentum j and therefore it is identically null if the sum of the 
spins Si + S2 is not an integer. 

3.1.3. Spectrum of distance operators 

At this stage, one can already ask the question of the quantization of distance opera- 
tors pip and their action on the generalized kinematical Hilbert space H-qCE, V). In 
fact, definitions of geometric operators have already been proposed in the framework 
of pure loop quantum gravity: these operators are self-adjoint operators acting on 
the kinematical Hilbert space and their spectrum have been computed [3 El 120] ■ 
It is claimed that the nature of these spectra gives the answer of the fundamental 
question of the space-time geometrical structure at the Planck scale. But this phys- 
ical interpretation is subject to polemics mainly because these geometric operators 
are not observables and therefore are not trivially related to a physical process. 

This section aims at clarifying this point, in the three-dimensional framework, by 
constructing classical distance observables labeled by a path 7 linking two particles 
()31|): we will identify the path and the associated distance observable in the sequel. 
Given two such observables, their Dirac bracket is non-trivial only if their intersection 
is not empty: if they share one particle, then the Dirac bracket is quadratic and 
involves a solution of CDYB equation [Hj; otherwise the Dirac bracket is computed 
using the Chern-Simons symplectic structure. 

In this article, we are exploring a quantization scheme based on loop quantum 
gravity techniques. However, the quantization is incomplete in the sense that we do 
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Figure 4: Distance observables and generalized spin-network states. The distance operator 
between pi and p2 is represented by a thick hne; holonomies defining generahzed spin- 
network states are represented by thin hnes. 

not quantize geometrical degrees of freedom living on the boundaries of particles. 
Therefore we cannot compute the action of a distance operator between two particles 
on a kinematical state involving one of these particles (left picture |3I); this case will 
be studied in [17]. It is nevertheless possible to make the computation in any other 
cases (right picture |^ which corresponds to what it is done in pure loop quantum 
gravity. 

Thus, let us consider a classical (square) distance operator characterized by 
a set of any two particles {pi, pj} and an oriented path 7 linking them. One can see 
that there is no ambiguity to quantize this observable and one obtains a quantum 
operator V'^{i,j) acting on 7io(S,'P). Let us consider a gauge invariant generalized 
spin-network state \E'5 G 7^0(^,7^) associated to a graph F such that there is one 
and only one edge C ^ ^^r which intersects the curve 7 on only one point x ^ ii,i2- 
It is clear that the action of T>'^^{i,j) on \E' reduces to the following action^: 

P^(^,j) > ^ = g^[e,A]tg^[e,A] > ^. (56) 

Let us recall that g^[e,yl] is the translational part of the Chern-Simons holonomy 
along 7 and its quantization is defined without any ambiguity by: 

q4e,A]'' = h^<^[A]^ie[ix)dx'' and q4e,A]l = h-^^lAfk eai{x)dx'' . (57) 
J 7 J 7 

Therefore, the result of the action ()56|) is an immediate consequence of the formula 
and reads: 

P?(^,J) > ^5 = s{x, C, 7)' ^c(^c + 1) ^ • (58) 

^Note that if the operator (hj) is clearly a Dirac observable, it is not the case for the operator 
on the r.h.s. of the equation (|56() . Therefore, this equation is not covariant and one could ask the 
question of its meaning. In fact, one has to understand this equation as a gauge fixing of the 
(square) distance operator and the r.h.s. of equation (|56|l corresponds to the expression of 
in a particular gauge. One can make such a gauge fixing in two different equivalent ways: one can 
make use of the translational symmetry to fix to zero the positions of the particles qi and qj or one 
can make use of the rotational {SU (2))-symmetry to fix the holonomy of the connection along 7 to 
a value such that any gauge transformation (translations in that case) is trivial. It is important to 
notice that, due to the relative configuration between the state and the distance operator, these 
gauge fixing do not affect the state itself. It would have been different if the state had involved 
the particle degrees of freedom in the operator. Finally, we see that the result of the action (|56|l is 
gauge invariant, thus well-defined and does not depend on the choice of the gauge fixing. 
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Figure 5: Action of the distance operator on a spin-network state. The action on the left 
is weakly equal to the action on the right. By weakly, we mean equal on the constraints 
surface, i.e. one can set any holonomy (with free ends) of the connection to be trivial on 
any contractible open set of S. 

This result is identical to the usual one obtained in pure loop quantum gravity even 
if our expression of the distance operator is a priori completely different than the 
usual one. Nevertheless, contrary to the usual (non-physical) length operator used in 
loop quantum gravity, this operator is well defined and one does not need to consider 
a regularization to compute its action. As a result, an important difference occurs 
when one considers the distance operator (j^UI) acting on states which intersect the 
path 7 several times (see left picture of Figure EI). In particular, one does not obtain 
anymore that the spin-network state is an eigenstate of the distance operator whose 
eigenvalue is a sum of contributions of different edges intersected. The action seems 
a priori more complicated and not necessary diagonal. However, as the distance 
operator is a Dirac observable, it should be diffeomorphism invariant and therefore 
should depend only on the homotopy class of the path 7. Thus, one could expect 
to find configurations where the computation of its action on a state is simpler 
and might be diagonal (see Figure Ej). To see this property, let us consider the 
example of the Figure El It is immediate to show that the equality between the 
two actions holds when one assumes that the holonomy of the connection along 
the part of the path 7 between xi and X2 is trivial, which is indeed the case on 
the constraints surface. Therefore, we have to make a gauge fixing which is state 
dependent: this point has to be related to the fact that the relation between gauge 
transformations and diffeomorphisms involves field-dependent gauge parameters. As 
a result, the action of the distance operator is invariant under local diffeomorphisms 
on-shell. Nevertheless, one has to be careful with global diffeomorphisms which 
does not generically leave the action of the distance operator invariant. This point 
is explained in the companion paper jTB]- Finally, one can really interpret colored 
edges of any spin-network intersecting 7 as fluxes of geometry and the invariance 
of the action of the distance operator under diffeomorphisms can be viewed as the 
conservation of the total geometry flux toward the path 7. 

Before ending this section, it is important to underline that we gave the spec- 
trum of distance observables acting on particular states. Given the distance operator 
between two particles, the contribution to its spectrum of generalized spin-network 
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states involving the degrees of freedom of these particles could be completely differ- 
ent. For instance, it seems clear that the action of 2) on the two-particle state 
will not be diagonal because of Heisenberg uncertainties. Indeed, such an action 
can be physically interpreted as the result of a measure of the distance between the 
two particles whereas the momenta of each particle are completely determined. The 
article [T7] aims to clarify these aspects. 

3.2. Spin-Foam amplitude 

This section is devoted to implement the remaining first class constraints (I35|) 
in order to define the Physical Hilbert space TiphysiX'^V) of the coupled system. 
The idea consists of generalizing the construction presented in ^B] by introducing 
a projection operator on physical states in the presence of particles. Then, we 
show that the physical scalar product admits a spin-foam state sum representation 
which defines a Ponzano-Regge model in the presence of massive spinning particles. 
Finally, we discuss the possibility to extend this construction to have a description 
of a background independent quantum field theory. 

3.2.1. Physical Hilbert space: projector and physical scalar product 

Given any function v G C°°(ZI,M'^), the functional $c(^) imposes constraints on the 
Dirac classical phase space of the coupled system. The constraints are first class 
if and only if the function v is constant on each boundaries otherwise they are 
second class 0. Thus, it is clear that the second class part of ^ciy) imposes that 
the component A^. of the spin-connection along each circle £j is constant and this 
condition has been already implicitly considered at the kinematical level. On the 
other hand, the first class part of the functional ^cl"^) imposes the condition that 
the spin-connection is flat everywhere on E and fixes its components A^p^ on each 
circle £j to the value: 

Vx e £i , A^X^) = A^X^i) = -^niiAiJoA-^. (59) 

At the quantum level, imposing these constraints on the generalized kinemat- 
ical Hilbert space should schematically identify spin-network states associated to 
graphs which are related by a diffeomorphism (connected to the identity). This 
is exactly what happens in the pure gravitational case ^Hl- Indeed, we construct 
an orthonormal basis B = of the physical Hilbert space TCphysiX') labeled by 

the so-called irreducible graph colored by a family J = (ji, ■ ■ ■ , jeg-e) of irreducible 
unitary representations of SU{2) {g > 1 being the genus"^ of the surface S). Any 
kinematical spin-network state labeled by a spin-network 5* = (F, (j-y)^, {iv)v) 
can be projected to the physical Hilbert space and its components in the basis B are 
given by: 

CIS = < ^5, >phys = < P*5, > , (60) 

^The case of the sphere and the torus are treated separately but do not introduce particular 
difficulties. 



24 









L 


L 




L l_ 









l_ 


(3 


(2) 

'€> 










r 










L 




L l_ 








1 











Figure 6: Cellular decomposition of the surface S with boundaries. For convenience rea- 
sons, we have drawn a regular square lattice. It induces natural simplicial decompositions 
of the boundaries associated to particles. 

where <,>phys is the physical scalar product defined from the kinematical scalar 
product <, > and the projection operator P on physical states acting on ?io(S) [TH] . 
To be more precise, the operator P is not well-defined in ?io(S) in the sense that it 
sends any kinematical state to an element of CyF(S), the algebraic dual of Cyl(S). 
Therefore, P is not strictly speaking a projector and the physical scalar product 
()6m) is in fact a bilinear form between Cyl(E) and CyF(S). Here, we have implicitly 
identified by duality the basis of Cyl(E) with its dual basis. 

From the definition (jHOI), we introduce an equivalence relation between kinemat- 
ical spin-network states as follows: 

~ ^5' if and only if V J , < ^5 - ^5', ^^J >phys = . (61) 

One says that the two states "^s and "^s' are physically equivalent. Given two 
spin-networks 5*1 and S2 such their associated graphs Fi and F2 are related by a 
diffeomorphism, one can show that "^Si ~ CH]- This property clarifies the above 
claiming about the fact that the constraints identify diffeomorphically equivalent 
spin- networks (as one could naturally expect). 

This section aims to generalize the previous construction to the coupled system. 
In order to define a projection operator on physical states, it will be convenient 
to consider, as in the pure gravitational cellular decomposition of the sur- 

face with boundary E such that each circle is contractible and does not contain 
any particle. The lattice intersects each boundary ii in z/j different points denoted 
{ipf'\ ■ ■ ■ , fi^'^) (figureini). As a result, it induces a decomposition of each boundary 
ii in z/j parts and we will denote by ef^^ the part going from ipf^^ to ipf^^^\ In 
the sequel, the set of edges and faces of the cellular decomposition will be respec- 
tively denoted by i5^(S) and iS^(E). To each element of e G 5^(S), we assign the 
5'f/(2)-group element he[A] and we will denote by Uf[A] = riee/^et^] ^ SU{2) the 
holonomy of the spin-connection around the face /. Using these notations, we define 
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the distributions P^. [y4,Aj] = HnLi -P^ where: 

(n+l) _ (n) 

Pi"^[A,A,] = 5(/i(„)[A]A,expK^^^-^Jo)Ar') • (62) 

' ATI 

It is clear that the distributions Pi. are a discretize version of the constraints ()59p . 
As a consequence, it seems natural to propose the following lattice version for the 
whole first class constraint $^7: 

TV 

P[A,(A,),] ^ U 6{Uf[A])l[PiAA,A,]. (63) 

There is no ambiguity in quantizing this operator for it involves only commuting 
variables. However, the obtained quantum operator P is clearly not well-defined in 
the generalized kinematical Hilbert space T-Co{Il,V) for the same reasons as in the 
pure gravitational case ^Hl- Therefore, P is not, strictly speaking, a projector but 
defines an application from ?io(S,P) to the space Cyl*(S,P), the algebraic dual of 
the space of generalized cylindrical functions ()52|1 . Nevertheless, we will identify in 
the sequel any element of a generalized spin-network basis in Cyl(S, V) to its dual 
generalized spin-network in Cyl*(S,P). Following the usual construction, we define 
the generalized physical scalar product between any two generalized spin-network 
states \l/5 and \I'5'' by the expression: 

< ^5, >phys = < P'^S, ^5' > • (64) 

Let us recall that <, > is the scalar product on the Hilbert space 7io(S, P). The ex- 
pression (j64p is defined for any cellular decomposition and therefore we can compute 
the scalar product between any two spin-network states based on any (continuous) 
graph immersed in the surface E. 

Two generalized spin-network states and "^s' are called physically equivalent 
(v]>_g, ~ vl/g,) if and only if they satisfy the following property: 

VF G noi^,V) , <^s-'^s',F>j,hys= 0. (65) 

To have a geometrical interpretation of this equivalence relation, let us consider two 
generalized spin-networks Si and 5*2 whose associated graphs Fi and F2 are related 
by a continuous diffeomorphism (p (connected to the identity). If the graphs do 
not touch the boundaries, then the spin-network states \E'si and "^32 are physically 
equivalent as in the pure gravitational case. However, we cannot generically conclude 
to the physical equivalence between the spin-network states when the graphs do 
touch the boundaries as we show in the following example. 

Let us consider a generalized spin-network Si whose graph touches (only) one 
of the circles ii in two different points (pi and ^i, the graph associated to the spin- 
network 5*2 intersects the same circle at the points and (figure [7j). The spin- 
network states "^Sa (q^ = 1) 2) associated to these graphs can be written as follows: 

^sAA,A,] = Plf") HKAA]):, hhdA])l, n\A,r:^ idjXndf ■ (66) 
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Figure 7: Example of two spin-networks whose graphs are related by a continuous difFeo- 
morphism (/>. The points ip\ and i9i are respectively sent to the points (^2 and '&2 by this 
diffeomorphism. The associated spin- network states ^ Sx and ^ $2 are physically equivalent 
to linear combinations of spin-network states associated to the graph on the bottom. 

In this expression, we have used the notation -R^jj," for the functions associated to 
the rest of the spin-networks drawn in (figure [7j). Let us recall that (j^-, i^J are the 
color associated to the boundary in that particular case, there is a unique choice 

3 k hi 

(up to a normalization) for the intertwiner Le- = leXJi ^jHi) '-V ® V ® V-^ C. 

Actually, the spin-network state Si is physically equivalent to a linear combi- 
nation of spin- network states s'^ whose spin- networks S2 differ from S2 only by the 
color associated to the boundary £j as it is shown in the following calculation: 

~ ^(/.,JA]A,e-^^o)^, ^(/.^JA]A,e-t^^")^,.,.(j,fc,j,jf 

~ ^(^JA]A,)^, HkAm^)l' e-^"^^^ idJ.k,JiX7 

In these lines of arguments, we have introduced the notations /i^a^i [^] (resp. /i,^2-,i [^]) 
for the holonomy of the connection along ^i between di and '§2 (resp. yji and (^2)'-, 
Ad = ^2 — ^1 and A99 = — V'l- This calculation makes use of the facts that 
t(j, fc, j^J is an intertwiner, h^^^^ [A\ and h^^^^ [A\ are physically equivalent to e™"' 2^"''' 
and e'"'^'^", and R^^b^ is physically equivalent to R^aif^ ■ The last line results from the 
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orthogonality relation ^ ., k,j'i_)L{j, k,j'f _) =7r(l)(g) 7r(l) and Cs^s'^ is a coefficient 
easy to compute. 

Of course, one can generalize the previous result to any spin-network state. In 
particular, any generalized spin-network can be physically extended as a sum of 
generalized spin-networks which touch the boundaries at their origin (figure Hj). For 
that reason, we will exclusively consider in the sequel latter types of generalized 
spin- networks. 

3.2.2. Ponzano-Regge model in presence of particles 

This section aims to picture the physical scalar product between any two spin net- 
work states s and s' in terms of spin-foam models. As in the pure gravitational 
case [TE], we will make use of this picture to construct, in a second step, an explicit 
basis of the physical Hilbert space. For convenience reasons, we will not give a gen- 
eral expression for the physical scalar product but we will illustrate the construction 
through typical examples which provide in fact a recipe valid for any situation. 

Before going to the coupled system, let us briefly review the main steps of the 
construction in the pure gravitational case [16]. Given two spin network states s 
and '^s' respectively defined on the graphs and F5/, we define the graph Vss' 
to be the minimal graph that contains F5 and F^/. To this graph is associated the 
notion of set of irreducible loops a^^ defined as the set of oriented boundaries of 
the corresponding graph. With these definitions, the physical scalar product reads: 

< ^5, >phys = < ^5, n -^(f^Tl^]) ^5' > • (67) 

The operator ^/^[A] denotes the holonomy of the connection around the irreducible 
loop 7 and acts on the states by multiplication. Then, we finish the calculation 
using the Ashtekar-Lewandowski measure. 

For the coupled system, the picture is similar but we use the expression 
for the projector. Thus, the expression of the physical scalar product between two 
generalized spin network states differs from the pure gravitational case only when the 
states involve particles degrees of freedom or when there exists an irreducible loop 
surrounding one or several particles. The following aims to show the construction 
of the physical scalar product in the latter cases. 

For that purpose, let us consider a state involving the degrees of freedom of one 
particle pi (figure |H1). The spin- network S associated to this state is such that its 
graph F5 touches the boundary ^i only at the origin and therefore separates the 
vicinity of the particle into different areas. Only one of these areas contains the 
particle and we will refer to it as the particle area. As a result, we define any spin- 
network state in such a way that the "virtual" particle is located in the particle area. 
Given two such spin-networks S and S", we can still define the set of irreducible loops 
a^^ in the same way as in the pure gravitational case. Nevertheless, we have to 
distinguish the subset af^ of irreducible loops surrounding particle areas from the 
subset a^^ of irreducible loops associated to empty areas. In that regularization. 
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Figure 8: Example of a generalized spin-network involving particles degrees of freedom. 
The vicinity of the particle is divided into three regions (a), (b) and (c) by the graph 
associated to this spin-network (on the right). The region (a) is called the particle area 
and the spin-network state (on the left) is defined such that its "virtual" particle is located 
in the particle area. 



the projector operator P reads as follows: 



P = U S{U,[A]) n 5(f/cJ^]A.e-^'«Ar')- (68) 

This expression of P is a trivial consequence of the definition (j63|) and Q denotes 
the irreducible loop surrounding the particle pi. This operator acts naturally by 
multiplication on any spin-network state. To compute its matrix elements, one 
proceeds as follows: first, one makes use of the SU{2) Peter-Weil decomposition 
formula of the Dirac distribution 6{x) = Ylj ^jXji^) iii terms of characters Xj(x) of 
SU{2) finite dimensional representations; then, one uses the generalized kinematical 
scalar product to finish the computation. The last point is performed using graphical 
techniques (see appendix of [TB]). 

To illustrate this construction, let us compute the physical scalar product be- 
tween the two states \l/5 and "^s' pictured in (figure 01) • The result gives the ampli- 
tude of the generalized Ponzano-Regge model in the presence of massive spinning 
particles. When one considers the normalization factors, the transition amplitude 
AI'^Sj ^5') between the states \l/5 and "^s' is given by the following formula: 

< P'^C \l/'c. > 

and can be pictured in terms of spin-foam model as shown in the picture (figure 
ITUI) . In that picture, one can see that the particle is associated to an edge but also 
to a face of the spin-foam. Indeed, the particle evolves along an edge but at the 
level of the vertex it chooses to follow one of the three emerging edges and therefore 
selects one of the three emerging faces (dual to edges). The selected face is called 
the particle face. 

In the case where there is no particle (m = and s = 0), the amplitude of the 
vertex reduces to the usual (6j)-symbol and we recover the well-known Ponzano- 
Regge model. In the case where the representations j and j' coloring the particle are 
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Figure 9: Physical scalar product between two states. The graphical rules are given in the 
appendix of ^S] and the box represents the group element e^"^" where m is the mass of the 
particle. The obtained result gives the amplitude of the Ponzano-Regge model associated 
to massive spinning particles. The presence of the dimensions Aj is just a matter of 
normalization. Note that the spin s of the particle is contained in the expression of the 
graph on the right (it labels the magnetic numbers of the open edges). 

trivial (particle without total angular momentum), the spin-network states do not 
explicitely involve the particle degrees of freedom but the definition of the projector 
P does: as a result, the amplitude of the vertex is still given by the (6j)-symbol but 
the amplitude of the particle face is modified. To be more precise, the face in the 
spin- foam colored by the representation r carries the weight Xr{^"^'^°) instead of the 
dimension A^: this face contains the information that there is a massive particle. 
This result is analogous to the one obtained in the covariant formulation of [TT] . 
Finally, one can define the model in the general case: the result depends on the 
choice of the intertwiners l and l', and the amplitude is given by a function labeled 
by 11 representations. 

At this stage, one can make use of the picture of the physical scalar product in 
terms of spin foam amplitudes to construct a basis of the physical Hilbert space when 
S is a Riemann surface of arbitrary genus g with an arbitrary number N of particles. 
In fact, the construction mimics the pure gravitational case [H^ and the result is 
illustrated in the picture (figure [TT|) for a surface of genus g > 1 with a number N 
of particles: we link all the particles together and we add the obtained chain to the 
graph associated to the pure gravitational case. One can easily see that any element 
of the basis is labeled by 6g—6+2N representations. Furthermore, one can show that 
any generalized spin-network state is physically equivalent to a linear combination 
of elements of this basis: in the spin-foam language, this property says that one can 
"evolve", by the action of the projector P, any spin-network state in T-Co(T,,V) to 
a linear combination of "simplest" spin network states based on graphs where the 
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Figure 10: Spin- foam picture of the transition amplitude between two generalized spin- 
network states. Note that, in that example, the particle face is the one colored by the 
representation r. 



maximal set of irreducible loops have been annihilated. Note that the choice of the 
basis is not canonical in the sense that there exist different parametrizations of the 
basis: for instance, one can glue the chain of particles to any edge (associated to 
pure gravitational degrees of freedom); or one can cut the chain in two parts and 
glue them on different edges... All these parametrizations give equivalent basis of 
the physical Hilbert space TiphysiX'^V). 

The case of the torus is similar: the basis is labeled by a graph obtained by 
adding to the pure gravitational basis graph the chain of particles. Any element 
of the basis is labeled by 2 -|- 2N representations. The case of the sphere is simpler: 
the basis is labeled by the chain of particles. When there is only one particle, the 
basis is trivial; when there are two or three particles, the basis is respectively labeled 
by 1 or 3 representations; when the number of particles is > 3, then the basis is 
labeled by 2N + 3 representations. 



4. Conclusion and perspectives 

In this article, we began by setting out the classical coupling between three dimen- 
sional Euclidean gravity and point particles. We have briefly presented the canonical 
analysis of the coupled system and we have expressed distance operators between 
any two particles as Dirac observables. Then, we have adapted loop quantum grav- 
ity techniques to quantize the coupled system: {i) we have generalized the notion of 
spin-network to provide the kinematical Hilbert space of the coupled system with 
an orthonormal basis; {ii) we have constructed a basis of the physical Hilbert space 
and the physical scalar product can be pictured in terms of spin-foam models. As 
a result, we propose a generalization of the Ponzano-Regge model which includes 
massive spinning point particles. Furthermore, the distance observables have been 
quantized and we have shown that their action on spin-network states reproduces 
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Figure 11: A basis of the physical Hilbert space for a surface of arbitrary genus g and a 
number N of particles on it. 

the well known result in three dimensional loop quantum gravity using non-physical 
distance operators. 

Even if we worked in the framework of Euclidean gravity without cosmological 
constant, most of our results can be easily adapted to the Lorentzian case and can 
include a cosmological constant. In all the cases, the classical dynamics of the cou- 
pled system is governed by the same action as the one we used but the gauge group 
is different and depends on the signature of the metric and the sign of the cosmo- 
logical constant [27J. Up to some technical subtleties one can define distance 
observables between two particles in a similar way. The quantization scheme is the 
same in the sense that one has to construct the kinematical Hilbert space and then 
the physical Hilbert space using spin-foam techniques. However, the analytic issues 
are slightly more involved in the Lorentzian case because of the non-compactness 
of gauge groups. As for the inclusion of a cosmological constant in the model, it 
should turn classical groups into quantum groups whose deformation parameter q 
explicitely depends on the cosmological constant (see for example). 

One of the main aspects of this article is certainly the quantization of distance 
observables and their evaluation on particular generalized spin-network states. Our 
distance operator is a well-defined operator on spin-network states and does not 
need any type of regularization, contrary the usual (non-physical) length operator 
used in loop quantum gravity. Its action on particular spin-network states have been 
computed and the result is slightly different than the usual one. All the same, the 
spectrum is still discrete. The algebra of distance operators has not been studied 
in details but certainly deserves a deeper understanding. In particular, it would 
be interested to compute the action of the distance operator between two particles 
on the associated two-particles state (figure 01): due to Heisenberg uncertainties, 
it is clear that this action will not be diagonal. Usual loop quantum gravity does 
not consider such situations and therefore it is not trivial to predict whether they 
contribute to a continuous or discrete spectrum. Whatever the result is, it will be 
certainly very interesting. 

A novel proposal providing a conceptual framework for background independent 
quantum field theories has been put forward in 0. The essential idea is the idealiza- 
tion of a measurement in a (background independent) quantum mechanical context 
by the specification of quantities on the boundary of a local space-time region. The 
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physical content in the boundary (e.g. its (space-time) geometry and matter field 
content) must be specified by a set (infinite in a theory with local excitations) of 
Dirac observables defined on the boundary Hilbert space. Because Dirac observables 
are very difficult to construct in four dimensions this proposal is presented at a for- 
mal level. Here we want to point out that a complete specification of the (physical) 
geometry of the boundary can be constructed in terms of the Dirac observables — 
introduced in section (2.2.2.) — that define physical distance between point particles. 
Of course there is no notion of physical time in this context as we are dealing with 




Figure 12: A system of point particles can be used as a net of Dirac observables to defined 
the physical geometry of the boundary of a local spacetime region. In the Lorentzian sector 
the distance operator can have eigenvalues corresponding to physical time. For suitable 
boundary states a notion of physical time elapsed between different particles leaving the 
space time region arises. 

Riemannian quantum gravity. In the Lorentzian context the distance operator con- 
tains time eigenvalues jTUI. There is no obstruction in principle to generalize our 
results to the Lorentzian sector. In that case preliminary results indicate that a 
notion of Dirac time elapsed between different events — defined in a diff-invariant 
way using the point particles — can be defined in analogy to the distance operators 
presented here. A representation of this idea is shown in Figure El A systematic 
study of this case will be presented elsewhere. 

Generalization to four dimensional gravity appears much more technically in- 
volved for several reasons. For instance, one cannot use point particles to define 
geometric operators in that case because they create a black hole. Then, it would 
be necessary to consider extended objects or matter field which complicate a lot the 
procedure. All the same, this point seems very important for it could put light on 
some open issues in loop quantum gravity as the meaning of the discreteness of area 
operators or the interpretation of the Immirzi-Barbero parameter. 

Finally, we would like to mention the possibility of making use of the previous 
results to describe an eventual self-gravitating quantum field theory. Let us recall 
that a consistent quantum field theory coupled to the gravitational field does not 
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still exist even in the three dimensional case. The main obstruction to define such a 
theory is that we cannot use standard quantum field theory methods: if one tries to 
quantize the coupling between the gravitational field and a scalar field for example 
using pcrturbative methods (developing around the fiat metric), it encounters that 
the theory is non-renormalizable. In fact, the philosophy of that method consists 
of quantizing first the field in a fiat spacetime and then considering the quantum 
fiuctuations of the gravitational field. Our result here provides a background inde- 
pendent quantization of a system of particles. As such it could be regarded as a 
first quantized version of particles and gravity coupled together providing the no- 
tion of n-particle Hilbert space 'H^^yg. There is no obstruction in putting together 
these Hilbert spaces in a second quantized version of the theory. Indeed, generalized 
spin-network states with N particles can be taken as an element of 

N 

As Fock space is the starting point for the definition of pcrturbative technique in 
standard quantum field theory, the previous Hilbert space could be of interest for 
the definition of generally covariant quantum field theories in three dimensions. To 
define the model, it is necessary to introduce a kind of interaction between particles 
that would mimic the self-interaction of a field. Therefore, this would modify the 
dynamics of the theory and the physical Hilbert space would be completely different 
than the one we have previously described. The implications of this viewpoint are 
currently under investigation. 

Avant tout... A nos Amours. 
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APPENDIX: Euclidean Group ISU{2) 

The Euclidean group G = ISU{2) is the isometry group of the three dimensional 
Euclidean space E'^. It is a six dimensional group and therefore any of its elements 
is parametrized an element z = {z^)i € M^; we will denote by X{z) the element of 
G parametrized by z. 

Its Lie algebra g is generated by infinitesimal generators of the three translations 
(Pi)i=o,i,2 and infinitesimal generators of the three rotations (Ji)i=o,i,2 which satisfy 
the following relations: 

[P„Pfc] = , [Pi,Jk] = e^k'Pi and [J,, J^] = e^k'Ji ■ (70) 

In these commutation relations, e*^' denotes the totally antisymmetric tensor and 
indices are lowered and raised by the Euclidean flat metric 6ik = diag(+l, +1, +1). 
We will generically denote by (^/)/ the six dimensional basis of g. 

The Lie algebra g is endowed with an non-degenerated invariant bilinear form 
(Killing form) <, >: g x g — > C defined by: 

< Ji, Pk> = S,k and < J,,Jk> = =< Pi,Pk> . (71) 

Up to a normalization factor, this non-degenerated Killing form is unique. Nev- 
ertheless, it is interesting to mention the existence of a degenerated Killing form 
(j ) : X ^ C defined as follows: 

iJ^.Jk) = and (P„Pfc) = = (P„Jfc). (72) 

Let us denote by t G g (S> g the Casimir tensor associated to the non-degenerated 
Killing form (fTTjl . i.e. tjj =< >; the components of its inverse are denoted t^"^. 
It is also useful to introduce the dual basis (^^)/ defined by the relation = t^'^^j. 

In this article, we identify the configuration space C of a relativistic particle with 
the group G and Euclidean transformations on the particle are given by the left 
multiplication Q: 

GxC ^ C , {g,X) ^ gX . (73) 

The left multiplication is completely determined by the application / : x — 
defined as follows: 

g{a)Xiz) = Xifia,z)), (74) 

where we have introduced the notation a = {(y^)i G and we have considered the 
element g{a) G G defined by g{a) = exp{—a^^i). In the article, we make use of the 
infinitesimal version of ()74j) which reads: 

By definition, the Euclidean group G is the semi-direct product of the group 
of three dimensional translations I by the group of three dimensional rotations 
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SU{2). In the vectorial representation, this property is expressed by the fact that 
any element X G G is represented by a 4 x 4 matrix as follows: 

In this representation, A is a SU{2) matrix evaluated in the spin 1-representation 
and q is an element of M^. We will denote by X = (A, q) this element and one can 
easily show that the group law is given by: 

(Ai,gi) ■ (A2,g2) = (A1A2 , Aigs + gi) • (77) 

As a result, the inverse element of X = (A, q) reads X^'^ = {A^^, — Ag). Given two 
elements X, F G G, we introduce the notation X1Y2 = X ® Y ^ G ^ G . To be more 
precise, we have Xi = X C?) 1 = XjEj 1 where {Ej)ij is the canonical basis of the 
space of 4 X 4 matrices. 

The group SU{2) is a sub-group of the Euclidean group. Unitary, irreducible 

representations of SU (2) are finite dimensional and labeled by a spin j G |N. Given 

j 

a half- integer j, we use the notation vr for the representation. The associated vector 
j 

space is denoted V, it is of dimension Aj = 2j + l and we introduce one orthonormal 

j 

basis {e^\m,n G [— j, These representations are also representations of the 
algebra su{2) and the basis is chosen such that it diagonalizes the element Jq G su{2). 
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